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Abstract 


Superconvergent  error  estimates  in  ^{H1)  and  t^H1)  norms  are 
derived  for  recovered  gradients  of  finite  difference  in  time/piecewise 
linear  Galerkin  approximations  in  space  for  linear  and  quasi-nonlinear 
parabolic  problems  in  two  space  dimensions.  The  analysis  extends 
previous  results  for  elliptic  problems  to  the  parabolic  context,  and 
covers  problems  in  regions  with  non-smooth  boundaries  under  certain 
assumptions  on  the  regularity  of  the  solutions. 
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1  Introduction 


1.1  Background 

This  paper  is  concerned  with  the  derivation  of  superconvergent  error  es¬ 
timates  for  recovered  gradients  of  finite  difference  in  time/piecewise  linear 
Galerkin  approximations  in  space  for  linear  and  nonlinear  parabolic  prob¬ 
lems  in  two  space  dimensions.  For  the  linear  and  quasi-nonlinear  problems 
the  analysis  covers  the  cases  of  variable  coefficients  and  of  regions  with  non¬ 
smooth  domains  under  certain  assumptions  on  the  regularity  of  the  solutions. 

Superconvergent  recovered  gradient  error  estimates  for  elliptic  problems 
having  solutions  with  low  regularity  due  to  Wheeler  and  Whiteman  [12]  and 
Goodsell  and  Whiteman  [4],  [5],  [13]— [15],  which  developed  the  work  of  Levine 
[7],  are  here  extended  to  parabolic  problems  with  variable  coefficients  and, 
using  numerical  integration  with  either  interpolation  of  coefficients  or  Gaus¬ 
sian  quadrature,  £2( H *)  results  are  produced.  The  analysis  has  been  further 
developed  to  produce  similarly  and  superconvergent  estimates 

for  quasi-nonlinear  parabolic  problems. 

The  application  of  gradient  recovery  for  time  independent  problems  was 
discussed  by  Krizek  and  Neittaanmaki  [6]  and  Wheeler  and  Whiteman  [12]. 
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Thomee  et  al.  [10]  have  considered  the  recovery  of  gradients  for  the  heat 
equation  in  domains  with  smooth  boundary,  and  in  this  context  have  proved 
£00(^00)  superconvergent  estimates.  The  results  here  involving  interpolation 
of  coefficients  extend  the  work  of  Douglas  and  Dupont  [1]  who  obtained 
optimal  convergence  results  in  for  nonlinear  parabolic  problems. 

1.2  Superconvergence  results  for  elliptic  problems 

The  results  for  the  parabolic  problems  of  later  sections  of  this  paper  are  based 
on  superconvergent  recovered  gradient  error  estimates  for  two  dimensional 
elliptic  problems.  These  are  reviewed  briefly  and  the  notation  is  defined. 

Let  ff  C  1R2  be  a  simply  connected  open  bounded  domain  with  polygonal 
boundary  dQ.  The  function  u(x)  satisfies 


-V  •  (fl(x)Vu(x))  =  /(x)  ,  xefl, 

(1.1) 

u  (x)  =  0  ,  x  6  dfl  , 

(1.2) 

where  a(x)  and  /(x)  are  continuous  and  their  second  derivatives  are  locally 
in  Too(fl)  and  a(x)  >  d0  >  0,  x  €  0.  A  weak  formulation  of  (1.1)  is  set  up 
by  multiplying  both  sides  by  v  6  and  integrating  to  obtain 

A(u,  v)  =  (/,  v)  v  veH^n),  (1.3) 
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where 

A(u,u)  =  f  aVu  •  Vu  ,  u,v  €  Hq(CI)  ,  (1.4) 

«/  n 

and 

(f,v)=[fv,  vzH‘(Sl).  (1.5) 

J  n 

For  convenience  and  ease  of  explanation  the  variable  coefficients  in  (1.1) 
have  been  taken  as  a  scalar  rather  than  as  a  positive  definite  matrix,  whilst 
homogeneous  Dirichlet  boundary  conditions  have  been  used  in  (1.2).  The 
matrix  of  coefficients  and  more  general  boundary  conditions  could  also  be 
treated. 

Following  [15],  it  is  assumed  that  for  the  region  0  of  problem  (1.1)- 
(1.2)  there  exist  subdomains  Do,  fl2,  such  that  flo  CC  fii  CC  fl2  C  0 
satisfying  the  following  conditions 

Condition  1.  ff0,  fii,  ^2  are  rectangular  and  are  partitioned  by  a  uniform 
isosceles  right  triangular  mesh  in  the  following  way:  each  subdomain  is  the 
union  of  a  finite  number  of  squares  each  of  which  is  subdivided  into  two 
triangles  T*,  by  the  diagonal  of  positive  slope. 

Condition  2.  fl2  is  remote  from  singularities  of  the  solutions  u  of  (1.3)  at 
corners  of  dfl,  so  that  u  €  H^(02),  and  the  regularity  of  /  is  assumed  to  be 
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such  as  to  guarantee  this. 

The  part  of  the  region  fl  —  CI2  is  covered  with  a  quasiuniform  trian¬ 
gular  mesh  which  is  compatible  with  the  regular  mesh  in  fl2-  Let  hk  de¬ 
note  the  diameter  of  the  escribing  circle  of  any  element  Tk  and  define  h  = 
maxejemen^s  hk-  A  finite  dimensional  subspace  Sh  C  Hq(Q,),  consisting  of 
continuous  piecewise  linear  functions,  is  defined  over  the  triangular  mesh 
partition  of  fl.  The  Galerkin  problem  approximating  (1.3)  is  that  of  finding 

Uh  €  Sh  9  A(uh,vk)  =  (f,vh)  V  vh  e  Sh  .  (1.6) 

In  (1.6)  it  has  been  assumed  that  all  integrations  have  been  performed  ex¬ 
actly.  If,  as  is  usual,  the  integrations  are  performed  numerically,  the  resulting 
problem  is  that  in  which 

K  €  Sh  9  A*(u*h,Vh)  =  (/,  vhy  V  G  Sh  ,  (1.7) 

where  A*  and  (•  ,  •)*  indicate  numerical  integration. 

We  wish  to  use  gradient  recovery  at  all  vertices  (nodes)  of  triangles  in  ST0 
and,  as  in  [4],  we  define  for  nodes  internal  to  flo  the  recovered  gradient  of 
Vh  G  Sh  to  be 

Hk{vh)=l-Y,{Vvh)n,  (1.8) 

D  i—i 
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where  T£,  i  =  1,2,  are  the  six  triangles  surrounding  the  node.  For 

nodes  on  the  boundary  dflo  we  define  suitable  specific  recovery  schemes, 
see  [13],  and  these  together  with  (1.7)  produce  recovered  gradients  of  Vh  at 
all  nodes  of  Oo-  To  these  we  now  fit  respectively  to  each  component  of  the 
recovered  gradient  a  piecewise  linear  function  (VvR)x  £  Sh  and  (Vu^)y  £  Sh. 
This  construction  ensures  that 

|Vu  —  Vu^lo.fio  <  C  (ju/  —  +  /i3|u|3,n0)  (T9) 

where  ui  £  Sh  is  the  piecewise  linear  interpolant  to  u.  The  following  result, 
which  uses  result  (1.9),  is  due  to  Whiteman  and  Goodsell  [15]. 

Theorem  1.1  Let  u(x)  £  Hq(LI)  fl  Hq(fl)  fl  i/3(02),  1  <  q  <  2,  be  the 
solution  of  the  weak  problem  (1.3)  and  umh  £  Sh  be  the  solution  of  the  Galerkin 
problem  with  numerical  integration,  (1.1).  IfVu*hR  is  the  recovered  gradient 
function  derived  using  (1.8)  then 

Vu  -  <  C(w)h^~"  (|u|,i0  +  (1.10) 

where  u  £  C°°(0)  fl  C^°((7)  is  a  cut-off  function  such  that  w(x)  =  1,  x  £  flo, 
u(x)  =  0,  x  £  (0  —  fix)  U  <9fh. 
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Remark  1.1  Result  (1.10)  has  assumed  the  regularity  u  €  Hq(Cl),  1  <  q  <  2, 
for  the  weak  solution,  thinking  of  lower  regularity  as  arising  from  the  data 
and  the  geometry  of  problems  (1.1)— (1.2).  It  could  be  that,  in  spite  of  the 
geometry  of  the  problem,  the  solution  is  in  fact  smoother  than  the  above  so 
that  u  €  H2(Q).  In  this  latter  case  the  estimate  (1.10)  becomes 


Vu  —  Vu 


.R 


o,n0 


<  C{u)h1+P  (| 


i2,n 


+  it 


3,n2 


). 


(i.ii) 


where  the  auxiliary  problem 


—V  •  ( aX7(f) )  =  if)  in  Q  , 
<f>  —  0  on  dVl , 

ip  6  is  such  that  <p  €  H1+p( fl),  0  <  p  <  1. 


(1.12) 


The  above  theorem,  remark,  and  recovery  of  gradients  will  now  be  applied 
for  the  space  derivatives  of  parabolic  problems. 
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2  Discrete  Time  Galerkin  Procedures  with 
Quadrature  for  Linear  Parabolic  Problems 

2.1  Weak  formulation,  notation  and  numerical  schemes 

In  this  section  we  consider  a  linear  parabolic  problem  with  nonconstant  co¬ 
efficients,  involving  the  elliptic  operator  of  Section  1.  For  this  problem  we 
derive  superconvergent  estimates  for  recovered  spacial  derivatives  of  discrete 
time  piecewise  linear  Galerkin  approximations.  These  estimates  will  be  ex¬ 
tended  to  more  general  linear  problems  and  in  Section  3  to  quasi-nonlinear 
parabolic  problems. 

The  function  u(x,  t)  satisfies  the  equation 

~  V  ■  (a(x,t)Vu(x,t))  =  f(x,t) ,  xeO,  teJ={0,T],  (2.1) 
together  with  the  boundary  condition 

u(x,f)  =  0,  x  £  <9f l,  (2.2) 

and  the  initial  condition 

u(x, 0)  =  'Uq(x)  ,  xGfi,  (2.3) 
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where,  as  in  Section  1,  tt  C  IR2  is  a  simply  connected  polygonal  domain  with 
boundary  dfl.  The  functions  a(x,  f),  /(x,  f),  (x,  t)  E  D  x  IR  are  assumed  to  be 
such  that  a(x,t)  €  Loo((0,T)-,  and  /(x,f)  is  positive  and  bounded. 

It  is  also  assumed  that  u(x,  t),  the  solution  of  (2. 1 )— (2.3)  is  such  that 

u(x,t)  E  n  L^H^n))  n  Hl(Hq(£l))  n  L2(H3(n2))  (2.4) 

where  l  <  q  <  2. 

For  convenience  we  have  assumed  homogeneous  Dirichlet  boundary  con¬ 
ditions.  General  Dirichlet,  Neumann  and  mixed  boundary  conditions  can  be 
treated  with  only  minor  modifications  in  the  arguments. 

In  approximating  problem  (2.1 )— (2.3)  we  use  a  continuous  piecewise  linear 
Galerkin  technique  in  the  space  dimensions,  whilst  the  time  discretisation  is 
treated  using  implicit  single  step  finite  difference  techniques  of  the  “0”-type, 
where  0  <  6  <  |.  These  time  discretisation  methods  are  chosen  because  they 
are  most  commonly  used  in  practical  computations  for  genuine  engineering 
problems. 

We  first  set  up  the  weak  form  of  (2.1)  by  multiplying  both  sides  of  the 
equation  by  a  test  function  u(x)  E  #o(D)  and  integrating  by  parts  so  that  we 
have  the  weak  form  of  (2.1 )— (2.3)  in  which,  for  any  t  E  J,  u(x,  t )  E  Hq(Q,)  x  J 
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satisfies 


(u,(x,t),v(x))  +  A(t;u(x,t),v(x))  =  (f(x,t),v(x))  ,  o(x)  € -f/j(S!)  , 

(2.5) 

where  ut  =  du/dt  and,  for  u(x)  €  Hq(Q), 


A(t;<f>,ip)  =  /  a(x,t)V(^(x))- V(^(x))dx, 

J  n 

(2.6) 

(/(  -  =  JQf(x,t)v(x)dx. 

(2.7) 

As  in  Section  1  we  again  define  Sh  C  Hq(Q)  to  be  the  space  of  piecewise 

linear  functions  on  fl.  The  continuous  time  Galerkin  problem  approximating 

(2.5)  is  that  in  which  Uk(x,  t)  €  Sh  x  J  satisfies 

(( uh)t( ■ ,  t),vh)  +  A(t\Uh,Vh)  =  (/( •  Vv*  €  Sh  , 

(2.8) 

where 

A(<;  u/!,  vh)  =  I  a(x  ,  l)V«/,(x  ,  t)  •  Vvhdx  , 

JQ 

(2.9) 

and 

(/(  •  ,  t),  vh)  =  /(x,  t)vh(x)dx  . 

(2.10) 

The  time  interval  [0,  T]  is  discretised  by  letting  At  =  T /N,  where  N  is  a 
positive  integer,  and  setting  tn  =  nAt,  n  =  0, 1, 2, . . . ,  N.  The  0-type  finite 
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difference  replacements  will  be  at  time  levels  tn'6  where 

tn’e  =  etn  +  (1  -  0)tn+l  ,  0  <0<\,  n  =  0,l,...,N -1  .  (2.11) 

Similarly  for  any  function  <f>(x,  t)  €  L2(ft)  x  L^J)  for  n  =  0, 1, . . . ,  N  —  1, 
we  define 

<t>n'd  =  6<j>n  +  (1  -  0)cf)n+1  ,  0  <  6  <  ]-  ,  (2.12) 

and  also  for  uf(x,  tn)  =  u"  the  backward  difference  replacement 

6,u  =^— = - Xt - '  (2'13) 

We  shall  adopt  the  following  standard  norm  notation.  If  w(x,t)  defined 
on  [0,  T]  x  U  is  sufficiently  smooth  and  1  <  p  <  oo,  then 

llu’ili,p((o1r),//)  =  !IF(0IIlp(o,T) 

where 

no  =  iMiww 

and  H  is  a  normed  linear  space  with  jj  j|^.  In  general  we  shall  take  H  to  be 
a  Sobolev  space  Hq{Vt)  for  some  positive  q. 

In  addition  we  define 

IML.oo.fi  —  IMlLooUo.rj.vy^fi))  (2-14) 


10 


iMloo.fi  -  IMMffO.TJ.LooCfi))  (2.15) 

where 

IHL.p.fi  =  ll‘ilwp»(fi) 

=  f  E  ll^-lll(fi) )  •  (2-16) 

\M<*  / 

For  p  =  2  we  have  that 

INI,  =  IML,2,fi  •  (2-17) 

For  the  corresponding  seminorms  we  use  the  notation  |  •  |  and  the  summation 
in  (2.16)  is  for  |a|  =  s. 

Lemma  2.1  Let  V  =  and  Q(-)  =  Q(x,t,u(x,t))  be  a  function  that 

is  bounded  above  and  below  by  positive  constants  ao  and  ai .  Further  assume 
that  Q  is  Lipschitz  continuous  as  a  function  of  u  and  that  |G|  and  |V-(5(u)G| 
with  G  =  ^  are  bounded  on  Cl  x  [0,T].  We  assume  the  elliptic  regularity 
given  by  (1.12).  Let  W  €  Sh  satisfy 

(Q(u)V(u-  W)  ,  Vv)  =  0,  v£Sh.  (2.18) 

If  u  and  belong  to  L2(0,  T;  Hq{Li)),  1  <  q  <  2,  then  u  =  u  —  W  satisfies 
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L2(0,X;L2(fi)) 


MIl2(0, T;L2(0))  + 


dv 

dt 


<  Kkq~1+P 

llU  lL2(0,T);//9(fi))  + 

du 

dt 

. 

L2(0,T,H<i(Q))\ 

Proof.  It  is  easy  to  show  that  there  exists  a  constant  C  such  that 


li,n 


<  Chq 


-l 


u 


?,n 


(2.19) 


Moreover,  if  a  6  V  is  such  that 


(Q(ii)Va,  Vu)  =  0  ,  v  6  Sh  ,  (2.20) 

then  it  follows  by  the  “Nitsche  lift”  (see  [9])  and  elliptic  regularity  that  there 
exists  a  constant  C**  such  that 

NU<C-4'IN|lfl  .  (2.21) 

To  obtain  estimates  for  we  employ  arguments  used  in  [11]  and  [2]. 

If  we  differentiate  (2.18)  with  respect  to  t  we  see  that 

(««)v§>)  .  (|£v,v.) 

=  (GVi/,Vv)  ,  v<=Sh.  (2.22) 
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One  can  easily  show  by  standard  arguments  that 


dv 

du 

dt 

<  Chq~l 

l 

dt 

?,n_ 

Let  /3  6  V  satisfy 

{Q(u)Vf3,  Vv)  =  (GVu,  Vv),  v€V  . 

Then  if  a  =  -  /3  we  see  from(2.21),  (2.22)  and  (2.24)  that 

IMIo.n  —  ||a||lin 

<  Chp+q~l 


r. .. 

du 

i — — 

IT 

■o 

b 

+ 

dt 

q,Q_ 

The  last  inequality  follows  by  noting  that 

Mhn<CMhn  ■ 

To  bound  ||/?||on  we  let  $  €  V  be  such  that 

(Q(u)VQ,Vv)  =  (P,v),  veV 


(2.23) 


(2.24) 


(2.25) 
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Then 


llfllln  =  (g(«)v*,v0) 

=  (GVi/,V$) 

=  (g(n)Vi/,g-1(u)Gv$) 

=  (g(n)vI/,v(g-1(u)G$-x)) 

-(g(n)Vi,,$v(g-1(u)G)),  xesh 

=  (g(«)vI/,v(g-1(u)G^-x)) 

+  (I/,v-g(u)*v(g-VK?)),  xes*. 

Thus, 

Iloilo, n  <  C’4’,M1flll*ll.+,fl  +  C  ll-ll,, nll*llx+p.n 

and  the  result  is  immediate. 

Returning  now  to  the  discretisation  of  the  parabolic  problem  (2.5)  and 
using  the  replacement  (2.13)  in  (2.5)  with  v  =  Vh  €  Sh,  we  have  a  t  =  tn,e 
such  that 

0 stun ,  t*)  +  u*)  =  (/(x,  r-9),  u*)  +  (£"’9,  vh) ,  vhesh,  (2.26) 
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where 

An’e{un'e ,  vh)  =  J  a(x,  tn’6)Vu(x,  tn'e)  ■  Vvhdx  (2.27) 

and  En'3  is  the  remainder  arising  from  the  time  discretisation,  i.e. 

En’e  =  8tun  -  ut{-  ,tn'd)  . 

We  are  finally  able  to  define  the  approximating  problem  for  (2.1)-(2.3) 
as  that  in  which  the  solution  u£(x) 

.  Sk 

satisfies 

toiO*)  +  A'n-‘  («;”■*.  v»)  =  (/(■■«"•*),»»)■  V116S*,  (2.28) 

where  A*n,e(-  ,  •)  is  the  numerical  approximation  to  An,6l(-  ,  •)  produced  by  us¬ 
ing  numerical  integration,  as  discussed  in  Section  1,  either  by  using  a  quadra¬ 
ture  rule  or  by  interpolating  to  the  coefficients  a(x,t)-  We  similarly  define 
(/(•  itn’8)ivh)*  to  be  the  numerical  approximation  to  (/(•  ,tn'9)-,vh)- 

For  every  t  6  J  we  now  define  the  elliptic  projection  Uh(x,  t)  €  Sh  x  J  of 
«  by 

A  ( t ;  t*fc(-  ,  f),  ufc)  =  A  (f;  «(•  ,  f),  ufc)  V  Vh  £  Sh  .  (2.29) 
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Then  u'l  =  Uk(x,tn)  is,  from  (2.26),  given  by 


(s,& ;,»»)  + 

-  A“’s  ((«,(.  ,("•'),»»)  +  (/(•  ,<"■*).  »*)  +  (£”•»,  v*)  .  (2.30) 

For  convenience  we  now  define 

(2.3i) 

r  =  (2.32) 

Then  subtracting  (2.30)  from  (2.28)  we  obtain,  using  the  above  notation, 

(Strin,vh)  +  A-n'6(Vn’\vh) 

=  (/(•  ,tn’6),vhy  -  (/(•  ,tn-V*)  -  (JS^.t,*) 

-  (^B,t,&)  -  A^(u^,Vfc)  +  An-*  («fc(.  ,*n’*),^) 

=  7?V)  -  (£"’%,)  -  (8te,vh)  +  T?'\vh)  ,  (2.33) 

where 

Ti,s(vh)  =  /(•  ,  (”■*,  %)•-(/(-,  <”•*),  t>») 

T2”'s(^)  =  A“*»(fi*(-  ,*“■•).«  -  . 
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Setting  vh  =  rjn’6  in  (2.33)  we  obtain 


But 


( c  n  n.6\ 

( 

[OtV  >7  ) 

= 

r/n+l  -  T]n 


At 


,drjn  +  (l-% 


n+1 


1 


,n+l 


=  a?!'1-'’* 

>  £{{<!-»)- 


o,n 
1  -20 


,n+l 


o,n 


O  +  ^llrf.n 


so  that 


>  2^{|h”+I|o,n-ll’>”llo,t.} 


and  hence 


2A  t 


,n+l 


-  Il-)"ll5.n}  +  b"’.  i'-’) 


<  +  |(E"V',)|  +  ■  (2.34) 

For  the  term  IT™'8 (r)n,8)\  in  (2.34)  we  have  that 


|TT V’*)|  =  !(/(•  ,  tn'e), -  (/(•  , t"'). t^) 


+n,6\  „n,0) 


<  «*!/(■  ,f-)|M{|^loja  +  h-1|wl}  (2-35) 
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using  approximation  theory  and  Lemmas  3.7  and  3.8  of  [12]. 


In  considering  the  term  \(En,e ,rjn’e)\  in  (2.34),  we  have  that 


\(En’e,r)n’e) |  <  C  En 


o.n 


o,n 


<  C 


{ni*iu  +  lh"+‘U  ■ 


Using  Taylor  series  expansions,  it  can  easily  be  verified  that 


En,e 

o,n 

<  C(At)2  ||nu||£,2((tnttn+l).  , 

(2.36) 

En,e 

o  ,n 

3 

<  C'(Af)2  ||Uttt||£,2((tnttn+l).  L2(0))  > 

(2.37) 

Now  consider  the  term  T£’e(r]n'9)  in  (2.34).  We  have  that 


r2n'V'1’) 


=  a  I”1')  -  a"'-‘  («;■*,  -)”■') 

=  A  (("'»;  u»(-  , -  A(tn'*\ u^X’) 


+  /!((* 


0 .  ~n,0  n,0 


;  u 


1 


(2.38) 


Consider  the  second  pair  of  terms 


A(tn,e;ul’e ,rjn’d)  -  A*n,e(u1’e ,t]n,e) 

=  J  (a(x,  tn'9)  -  a*(x,  tn'9))  Vunh’9  ■  Vr}n'9dx  . 
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For  the  case  of  interpolation  where  a*(-  ,  •)  =  a/(-  ,  •)  this  gives 


a  /  j.n,8  ,»  *n,0 n,0  „n,0  \  |  /iL  2  ||  || 

A(t  ,nk  ,T]  )—A  (uh  ,T)  )|  <Ch  II a ll2,oo,« 


-  n,0 

71,0 

1,0 

V 

l.n 


(2.39) 


For  the  quadrature  case 


£  [  (<■(',  '“•*)  -  a(Gt,  (”•*))  VfiJ'*  •  V?*dx 

I.  JT'j, 


<  C h2  II  a  I 


2,oo,n 


~n,6 

71,0 

uh 

1,0 

V 

i,n 


(2.40) 


Now 


~n,S 

Ul’ 


1,0 


< 


u 


n,8  n,B 
—  U 


1,0 


+ 


u 


n,0 


i,n 


<  c 


u 


n,0 


ii.a 


+ 


£71,0 


l,n 


(2.41) 


Thus  from  (2.38)-(2.41)  for  t  €  (<n,tn+1)  we  have 
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For  the  first  two  terms  in  the  right  hand  side  of  (2.38)  we  have 


A(r*  ;«*(• . 

=  I  o(x,  ,  i"-«)  -  <•') . 

J  n 

=  jf  a(x,r-»)v  («»(. ,  ("■*)  -  «s;  -  (1  -  «K+1)  •  v^-Ux 


<  llal 


<  < 


00,17 


r,n+J 


,n>9 


i,n 


li, n 


loo.n 


loo.n 


,n,e 


-TlyO 


10  X  CAt*  |(tt*)t|L2(((nitn+l);  tfl(ft))  ’  ^  ^  <  2 

IQ  X  K^)«lL2((t",t"+l);  Z/1  (n))  >  ^  =  2  • 


(2.43) 


Estimates  (2.42)  and  (2.43)  together  provide  a  bound  for  IT?'9 {r]n'e)\. 
Consider  next  the  term  (6t£n ,T]n'e)  in  (2.34);  by  Cauchy-Schwarz 


Now  =  §i{uh  ~  u")  and  fr°m  Lemma  2.1 
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<  A'-,+'C’(||<>||,», n,||a, |L,„  X 

MlL*((i",tn+1);tf«(n))  +  llU*ll/,2((X'Vn+1);tf,(«))) 


io,n 


+ 


n+l 


0,Q 


(2.44) 


From  (2.34)  and  the  estimates  (2.35)-(2.44)  we  find  that 
1 


n  4*1 


2A  t 

<  \ch2\\f(-,r-(‘) 


2,n 


+  h’-'*rc  (lhlL,n  ll»IL.n) 


u  II  L2((t",e*+i );//«(«))  +  Hu*llL2((t’>,t"+>);m?(n))^ 
^+1u| 


dtP+' 


+  CAT3'2 

II  yifWtl  I, 

IZ,2((eV"+>);<T,2(fi)) 

+  C  ||a|L,n  (Aff-i 


o,n 


+ 


,n  +  l 


o,n 


+ 


d0uh 

} 

dt0 

L2((«n+tT,+  l);Hl(n))  J 

\rin'6  |i,n  (2.45) 


where  /?=lifO<0<|  and  /?  =  2  if  6  = 

Let  do  =  inf(o,r)xn  a*(x, 1)  and  fli  =  sup(0 j)xn  a*(x,  t).  Then 


2 

o  ,n 


(2.46) 
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Thus,  using  (2.46)  and  the  inequality  ab  <  ^-a2  +  |62  with  e  =  |  for  terms 


involving  || 77" ||0  ^  and  !|i?n+1|lon  and  e  =  a0  for  terms  involving  r)n,e 


1,0 


we 


obtain  from  (2.45)  that 


2 

0,0 


<  Ch‘  ||/(.  +  C  (W|„  n  ,  ||a,IU,n) 


X  <  u 


lll(L2(t’V"+1);  H*(Q))  +  llU«ll(L2(t’\t"+1);//«(0))} 


+  C(Af)2/3_1 


dp+] 


u 


dt 0+1 


Z/2((<n,tn+1 );  L2(0)) 


+  C(a0)  ||«Hl«(0,T);  ^2,(0))  I  ^ 


U 


n,6 


1,0 


+  At20-1 


dpu 


dte 


Oo 

+  2 


n,8 


1,0 


+  C 


(llxlo.n 


+ 


Ls((t»  ,*"+»);  Hl{ 0)) 

2  ) 

0,0  /  ’ 


n+1 


(2.47) 


where  in  (2.47)  we  used  the  fact  that 


u 


n,6 


1  n  <  C(a0,a j) 


u 


n,6 


1,0 


and 


dpuh 


dt? 


< 


0,0 


C  (a0,  llallw^to)) 


dp 


u 


dtp 


1,0 


Multiplying  (2.47)  by  2At  and  summing  over  n  for  n  =  0, 1, 2  . . . ,  N  —  1,  we 
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obtain  the  estimate 


N 

2 

0 

n 

0,0 

V 

N- 1 


*|-  2a0  b  At 
0,0  “ 


n=0 


,n,( 


1,0 


'iV-l 


2 

2,0 


<  C'1/i2('J_1+p)  t  ]P  At  ||/(-  ,  t”,e)| 

\  n=0 

+  c,(IHU5IKIIoo,o)II  ull(//i(o,T);//9(n)) 


+  IMIz,°°(0,T);  W^O))  llUll(L2(0,T);//l(O»  f  + 


+  C2(At) 


2/3 


d0+] 


u 


dt0+l 


L2(0,T);  L2(0) 


+  «l 


(L~(0,T);1V|,(0)) 


dtp 


(L2(0,T);  , 


+  &  { E  Ai  Itolio  +  A(  lhA'lC„} 

l  n=l  ’  J 


(2.48) 


which  can  be  written  as 


(1-CAt)  rjN 


N—l 


0,0 


+  do  ^  At 


n=l 


,n,0 


1,0 

/V-l 


<  C!/i2(9-1+p)  +  C2  (At)20  +  C3  £  At  ||i?n||oin  (2.49) 


n=l 


where  the  constants  C\,C-2,Cz  depend  on  various  derivatives  of  f,A  and  u 
as  in  (2.37). 

Application  of  Gronwall’s  Lemma  (see  [2]  and  [3])  to  (2.49)  yields  the 
following  result  which  we  state  as  a  theorem: 
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Theorem  2.1  Letu(x,t)  G  L00(Hq(Q.)C\HQ(Cl))r\H1(Hq(Ll))  fl  L2{H2’ (Vt2)) , 

l  <  q  <  2  be  the  solution  of  the  weak  form  (2.5)  of  the  parabolic  problem 

(2.1)-(2.3)  defined  in  fl  x  J.  Let  Sh  G  #o(fl)  denote  the  spaces  of  piecewise 

linear  functions  defined  on  the  triangular  partition  of  fl  as  in  Section  1  and 

*n 

let  Uh(x,t)  G  Sh  be  the  elliptic  projection  ofu(x,t).  If  uk(x)  is  the  solution 
of  (2.28)  at  tn ,  where  J  =  (0 ,T),  T  =  NAt  and  tn  =  nAt,  then 


*N  ~N 
uh  ~uh 


2  N-l 

+  o,q  ^  At 

O.fi  n=0 


(1  —  At) 

<  Clh2^-l+p)  +  C2(At)213 


<  n,6 


Uh  ~uh 


n,9 


i,n 

(2.50) 


where  C\  and  C2  are  constants ,  l<q<2,/3  =  l  when  0  <  9  <  |  and  (3  =  2 
when  9  =  Here  0  <  p  <  1  is  a  parameter  arising  as  in  (1.12). 


At  any  time  level  tn,  n  =  1 ,  we  now  apply  the  recovery  proce- 

*n 

dure  described  in  Section  1  to  uh  to  produce  the  recovered  gradient  function 
V  uh  .  Thus,  combining  the  inequality  (1.10)  of  Theorem  1.1  with  the  re¬ 
sult  of  Theorem  2.1,  we  now  have  the  main  (superconvergence)  result  of  this 
section. 

*U 

Theorem  2.2  If  u(x,t)  and  uh(x)  are  as  in  Theorem  2.1 ,  and  all  the  hy- 

*  n,R 

potheses  of  that  theorem  hold,  and  i/V  uh  is  the  recovered  gradient  function 
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at  the  time  level  n,  n  =  1, 2, . . . ,  N ,  then 


where  1 
where 


Vu  —  V  u 


R 

h 


2 

(L2(0,t*);L2(n0)) 


=  0{h2{q~l+p))  +  {At)20 


(2.51) 


<  q  <  2,  0  <  p  <  1,  and  (3=1  if  0<6<\  and  /3  =  2  if  9  = 


lll^7u,lll(i:2(0,tJV);L2(fio))  —  At\Vu(t  )| 


71=0 


2 

LHQo)  ■ 
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3  £2(0,  T)  Recovery  of  Gradients  for 


Nonlinear  Parabolic  Problems 

We  now  extend  the  previous  analysis  to  nonlinear  parabolic  problems  in 
which  the  function  u(x,t)  satisfies  the  equation 

Qu(x. 

— 7^ - V  •  (a  (x,  t,  u(x,  <))V«(x,  t))  =  f(x,t,u(x,t)) , 

xeO,  t  e  J  =  (o,T] ,  (3.1) 

together  with  the  boundary  condition 

u(x,t)  =  0,  x  €  dQ  ,  (3.2) 

and  the  initial  condition 

u(x,  0)  =  u0(x) ,  x  6  fi  .  (3.3) 

The  same  regularity  conditions  on  /  as  in  Section  2.1  are  taken,  and  we 
further  assume  that  /  and  a  are  Lipschitz  continuous  with  respect  to  u  with 
Lipschitz  constant  K ,  and  that  the  regularity  of  u(x,t)  is  as  in  (2.4). 

The  weak  formulation  corresponding  to  (2.5)  is,  for  any  t  €  J,  to  find 
u(x,t)  €  Hq(Q)  x  J  such  that 
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(ut(x,i),u(x))  +  A(t;u(x,t),v(x)) 


=  /(x,t,u(x,t),u(x))  Vw(x)  €  fTo(fi)  ,  (3.4) 

where 

A(t;<f>,ij>)  =  f  a  (x, «(x, <))  V^(x)  •  V^(x)dx  (3.5) 

Jn 

and 

(/(•,f,u(*,<)),v)  =  J^f{x,t,u(x,t))v{x)dx.  (3.6) 

Equation  (2.26)  now  becomes 

(6tUn,Vh)  +  A(tn'e-,u(;tn'e),Vh)  =  (f(x,tn’*,u(x,tn’*),VA) 

+  (En’\vh)  VvheSh  ,  (3.7) 


whilst  the  Galerkin  problem  with  quadrature,  corresponding  to  (2.28),  is 

(3.8) 

where  the  A*n,d  is  now  defined  as 


A'n'e{<t>n'\  vh)  =  [  o*(x,  <n’*,  Ex<j>n'9)V<j>n'e  ■  Vvhdx  .  (3.9) 

J  f) 
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Here  the  use  of  the  *  in  a *  and  (/,  •)*  has  the  same  definition  as  in  Section  2 
and 


E‘<pn'“  =  ^  ,  0  <  0  <  j 

=  ,  «  =  i.  (3.io) 

In  the  analysis  below  we  denote  by  //  the  linear  interpolate  in  Sh  of  /. 
Equation  (2.30)  correspondingly  becomes 

(6lul,vk)  +  A"''s(i"K-,,vk) 

+  (/(-,  r*,U(x,  (“'»)), »»)  +(£“•*,  vh).  (3.11) 

Modifications  to  extend  the  result  of  Theorem  2.1  to  the  current  nonlinear 
problem  involve  only  the  reestimation  of  T”’  and  T^’6  for  the  present  case. 
In  order  to  do  this  we  proceed  as  follows.  Following  (2.24)  we  now  have  that 

7?-V)  =  (/(x,  tn’e,  Exu™’e),  Vh)* 

-  (f(x,  u(x,  (”■*)),  (3.12) 

Let  {zj(a;)}^1  be  a  nodal  basis  for  Sh  corresponding  to  the  set  of  nodes  Xi 
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of  the  partition  of  Q,  and  define 


M 


fl(x,  t ,  tt(x))  =  u(xj))zs  • 

i=i 

Set  \4  =  T]n'e .  We  first  consider  the  case  of  interpolation  of  coefficients.  Thus 


T?’e(r]n’e) 


X, E*«rV<))  -  /(Xy,  (“•»,  ufo,  ("■•)))  Z,,  ,"•< 

+  (/,(x,  <"•*,  u(x,  ("•*))  -  /(X,  U(X, 


< 


E  g£(x>,*“Vi)  (fi'CV;)  -  U(x„r’#)2i, 

3 

+  Ch 2 1|/(("'»)|2  n  ||,-« 


0,0 


where  ^  lies  between  Exu™'e(x.j)  and  u(xj,tn'9). 
Now,  assuming  that  |<9//<9u|  <  K,  we  have  that 


E  ^(x,-,  -  u(xj, 


< 


KY.\E’<'‘(x,)-E-uY(x,)\z, 

J 

+  '£\Exuf(x,)-u,(xJ,r'e)\zi  . 


(3.13) 


Thus, 
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where  /?=1,O<0<|  and  (3  =  2,  6  = 

Here  we  have  used  the  equivalence  of  norms  for  finite  dimensional  spaces, 
i.e. 

We  simplify  the  bound  (3.14)  as  follows: 
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rr,l’(i)n',)|  <  c(i<)  |  ||ij"|i5i0  + 

+  iinofl  +  |r-,|C(1  +  ii(«/-“riio.i! 

+  |(U/-ur'i;n  +  (Ai)“pP)|[n 

+  if'6  2  )+Ch4  f(tn'e)  n 

1  o,n  j  J  v  '  2,0 

<  C(K)  |  +  I’)""1  £q 

+  hn,-^>  £  yf  +(At)»  §^(iV  ) 

(=n-l  ’  Ul  1,0  J 

+  C/.4jj/((“'»)||2n  +  C||,”+1|[!1  (3.15) 

where  tn  is  a  point  between  tn~l  and  tn,e . 

Now  using  Gauss  quadrature  we  have 


i3 

3 

Qa 

II  ll 

E/  [/(Gt,l"'«,£*ur"(Gt))~  /(x,i"'V(x,  ("•')) 

k  JTk 

E  f  [S,Y  +  sji*  +  sg  +  sj;']  -i”'^ 

(3.16) 

where 

nn,0 

^lk 

=  f(G„,  ("•«,  £'«r''(Gt))  -  /(<?*,<”•*,  E‘uf(G„)) 

(3.17) 

QTijO 

k 

=  f(Gk,C\E‘uY(Gk))  -  SW\E’uY(x)) 

(3.18) 

QTlyO 

^3  k 

=  /(x,  ("■*,  £*«?•(!))  -  f(x,  «,(*,  <"■*)) 

(3.19) 
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SnJ  =  f{x,tn\Ul{xX'6))- f{x,tn'6,u{x,tn'e))  (3.20) 

Using  the  Lipschitz  continuity  of  /  and  the  definition  of  Ex ,  we  have  that 

L  /  (sa*  +  s4Y)-)"’#«ix 

^  .m 


<  c 

(At)20-1 

d0u 

dtP 

2 

+  K  llu  —  ui  Iloo((o,t),l2(o)) 

. 

L2((tn,tn+'),L2(n)) 

+  C 


o,n 


(3.21) 


Now  for  simplicity  we  consider  first  S^'k  . 


E  L  Stf^dx  =  J2  L  sXV'  -  Vn'e(Gk)) 

K  JTfc  L. 


Tk 


1  ,Tk 


<  £/>!t||/(x,r^.EX'V))|1J|^ 

£  »o  +  C/>4  ||/llh,  (3.22) 

In  obtaining  (3.22)  we  have  used  the  fact  that  for  any  continuous  F(x) 

/  (F(x)  —  F(Gk))  ■  Idx  =  0  , 

JTk 

thus, 

[  (F(x)  —  F(Gk))0(x)dx  =  [  (F(x)  —  F(Gk)) (6 (x)  —  const) dx  . 

JTk  JTk 


Now  turning  to  5”/,  we  have  that 
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E  L  stf d* 

kjT> 

=  E{/r  [f(GkX'e,Exuhn’e(Gk))  -  f{Gk,tn\ExuY{Gk)) 


-  /(x,  E*uhn'\x))  +  /(x,  E*Ul(x)) 


^n,$ 


+  J  [/(x,^n•^^u^fl(x))-/(X,<n•^^u/(x))]^dx  (3.23) 

By  applying  the  mean  value  theorem  and  Gauss  quadrature  estimates  to 
(3.23)  we  have  that 


E  L  s^dx  <  EMkiL  | 

Ic  JTk  t 


o.r* 


+ 


i  r r  771X  *ri, Q  t-ix  71,# 

+  A  Euh  —  E  Uj 


„n,0 

o,n 

»7 

o,n 


-  Elk 

k 


n,6 


1  ,T* 


5/ 


rpx  tx,0  t?x  *n,0 

E  uj’  -  E  uh’ 


1  iv'  rr>x  *n.$  t^x  ti, 

4-  A  E  uh  —  E  U[ 


Looqo.r^LoofQ)) 

0,Q 


0,T* 


n,0 

o,n 

V 

Thus 


E^  5in^x  - 


.,n,0 


i,n 


+  C 


+  c,in>?niiin  + 


.n  — 1 


+  c  E  -  u/) 


f=n— 1 


2 

o,n 

2 

o,n 

2 

o,n 


(3.24) 
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Combining  (3.21),  (3.22)  and(3.24)  we  have  an  estimate  for  Ti'e(r)n’9).  The 

only  terms  which  did  not  appear  (for  the  linear  problem)  in  (2.45)  are  the 

2 

il7/"-1 1 1  o  n  and  (^h  -  UlY  ,  t  —  n  —  1,  n,  in  (3.24).  These  of  course  arise 
from  the  extrapolation  which  was  not  needed  in  the  linear  case.  This  has  no 
effect  on  the  final  estimate,  as  in  (2.45),  because  of  the  application  of  the 
Gronwall  inequality  and  approximation  theory. 

We  now  seek  to  estimate  the  term  T£'9(T)n,<>),  as  in  (2.27),  for  the  nonlinear 
context.  The  term  T%'6{i)n'6)  is  defined  as 

T^e{r]n’e)  =  -A-n’\unh'9,Vn'e) 

+  />((”•*;  i»(-,  !"•»), 


In  a  similar  manner  to  the  treatment  of  for  the  case  of  the 

interpolation  of  coefficients  we  have  that 


x2n’V,<?)  = 


=  -/  E«(xi,^,^r<,(xJ))zJvur-v^dx 
Ju  j 

+  /  a(x,r'e,u(x,^))Vu,(-,r’e)- v^dx 

J  n 

=  /  +  G(xj,^,X^(xi))zjV^  •  Vrf’*dx 

Ju  j 

+  [  V  -«(xj,  (”•«,  fi'uf'fxj)  +  a(x„  (”■*,  U,(x„  tn,e))zjVu^'e  •  V,"’Vx 

lai 

+  Ja  fe -a(xJ,i"’‘’,u,(xJ,  (“•*))*,  +  a(x,("'e.u(x, «”•“)) j  •  V,”'"* 
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+  J  a(x,tn'd,u(x,tn'e)(-Vunh'e  +  Vufc(-,  <"**))•  Vijn'8dx. 

The  proof  proceeds  similarly  to  the  treatment  of  the  T^'6{ijn'6 )  terms.  In 
particular 


T?\ rtn'e) 


<  I< 


+ 


E  EX<n'\x3)  -  EV/iXj) 


o,n 


o,n 


+  C/.2||a||j  n  (<”'«)  |  ||Vfit||„i0  ||V,”'‘ 
I  d0u 


+  C(At)0~ 2 


dtt> 


o,n 


o,n 


so  that 


<  WHVfilUoJx 


lignin  +  IK1 


+  (At)2/3_1 


dt* 


o,n 

2 


2 

2, ft 


+  /j2(g-l+p)  ^ 

i=n— 1 

+  a4  Nlf0  («"'•) 


L2((t",t"+1);H‘(ft)) 


+ 


do 


71,0 


0,ft 


(3.25) 


35 


We  remark  that  for  q  =  2 


liviuiu 

£  llv“IL.n  +  h~'  !lv(“/  °°(0,T);L2(n) 

+  l|V(ti-ti/)||00i0 

=  0(1) 


since 


l|V(u/  ^/i)IIl^({0,T);L2(Q))  —  hq 


(3.26) 


In  the  case  of  Gauss  quadrature  the  treatment  of  T£’e(rjn'9)  simplifies 
since  Vu£,<?  •  S7r]n'9  is  constant  over  T*.  More  precisely,  we  have 

T^V’*)  =  Y,(-  f  a(GkJn'd,E*ur’e(Gk))Vunh'e -V^dx 

TV  Jt> 

+  [  a{x,tn'9,u{x,tn'9))Vunh'9  ■  Vvn’9dx) 

JTk  / 

=  E  (-  /n  «(x.  <"•*,  £'«r,#(x))  vu;'#  •  v,"'Vx 

+  /  a(x,f's,«(x,<n’e))V«^-  Vrin'edx 
J  Tt  ' 


Vul 


oo,f2 


Vr/ 


n,0 


0,fl 


.(3.27) 


The  bounds  for  T™'8{r]n’B)  and  T2”’e(?/n’9)  include  no  new  terms  over  those 


for  the  linear  case  except  the  ||^n  1||0>n  and  YJi=n-\  (“fc  —  ui ) 


o,n 


so  that 
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an  estimate  of  the  type  (2.47)  can  be  obtained  and  this  leads  immediately  to 
the  theorem 


Theorem  3.1  Let  u(x,t)  be  the  solution  of  (3.  l)-(3.3)  and  have  the  same 
regularity  as  that  given  by  (2. 4).  We  further  assume  with  a(t)  =  a(t,u(t)) 
the  elliptic  regularity  assumption  for  q  —  2.  Let  u*hn(x)  be  defined  by  (3.8) 
with  =  u°h  or  uj(-,  t°).  Then 


N-l 


(1  -  At)  7]N  oq  +  Go  Y,  At 


n=l 


.n, 6 


<  Cxh2(q~l+v)  +  Ci(At) 


20 


2 

1,0 


where  (3  =  1  when  0  <  9  <  |  and  (3  =  2  when  9  = 

In  addition  the  following  theorem  results  immediately  from  Theorem  3.1, 
and  the  recovery  results  of  Section  2. 

Theorem  3.2  Let  u(x,t)  and  u*h(x,t)  satisfy  the  hypotheses  of  Theorem  3.1 
and  let  Vu*hnR  be  the  recovered  gradient  function  at  the  time  level  n,n  = 
1,2  ,...,N.  Then  (2.51)  holds  for  (Vu  —  \7u*hR)  in  the  case  of  the  nonlinear 
problem  (3.1)-(3.3). 
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4  £oo(0,T)  Recovery  of  Gradients  for 

Nonlinear  Parabolic  Problems 

All  the  estimates  derived  so  far  in  this  paper  have  been  in  terms  of  the  L2- 
norm  in  time.  From  a  practical  point  of  view  it  is  much  more  useful  to  have 
estimates  which  in  time  are  in  terms  of  the  4,0-norm.  We  now  derive  an  £0 0- 
norm  time  estimate  for  the  backward  difference  procedure  ( 0  =  0)  as  applied 
to  the  nonlinear  problem  (3.1)— (3.3).  Whilst  for  simplicity  we  consider  the 
case  0  =  0,  our  arguments  can  be  generalized  in  a  straightforward  fashion  to 
include  the  case  0  <  0  <  |. 

We  proceed  as  in  the  previous  sections  and  taking  0  =  0  subtract  equa¬ 
tion  (3.11)  involving  the  elliptic  projection  uk  from  equation  (3.8)  involving 
the  calculated  Galerkin  approximation  u*h.  For  ease  of  notation  in  allowing 
variability  of  the  third  variable  in  a( ■  ,  •  ,  •)  we  redefine  the  A*  and  A  to 
contain  three  arguments.  Thus 

A*n+1(z,  <f>,  ij>)  =  a*(x,  tn+1,z)V<j>  ■  V^cfx 

/  a/  (x,tn+1,z)V<f>  •  V^cfx 

=  -  (4.1) 

W  a(Gk,tn+\z(Gk))V<f>-V^dx 

.  k  jTk 


38 


depending  on  whether  interpolation  or  Gauss  quadrature  is  used,  and 


A(t-,  z,  <f>,  ip)  =  I  a(x,f,.z(x,t))  V<f>- V?/>dx  (4.2) 

J  fi 

so  that 

=  -  A-”+VrX+\»<.)  +  A-"*l(u\ul*\vk) 

-  («“,»»)-  Am*'[u',a?\»k)  +  A((“+,;u“+1,iJ+1,»*) 

+  (/(x,i-+',Urtl),t.k)-  -  (/(x,  (”+1,«n+I),»0  -  (£"'%*) .  (4.3) 

Setting  Vh  —  StT]n  in  (4.3),  multiplying  the  result  by  At  and  summing  the 
result  over  n,  n  =  0,1,2,...,  TV  —  1,  we  obtain 

£  A(  Kilo.n  +  E 

n=0  n=0 

N- 1 

=  X]  +  ^2"  +  ft) 

n=0 

-  £  At  [(*t£V«»?n)  +  (£"’*, <W')]  (4.4) 

n— 0 

where 

ft  =  -A*n+1«n,<+1,^n)  +  A*n+1(un,uj;+1,  W1)  ,  (4.5) 

ft  =  —  A*n+1(t*^,«S+1,^t7n)  +  A(tn+1;tin+1,«S+1,^n)  ,  (4.6) 
ft  =  (/(x,t"+1,ur+1),^)‘  -  (/(x,tn+1,un+1),^n)  .  (4.7) 
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Writing  the  time  difference  8tr]n  as  ( rjn+l  —  ijn)/At  and  taking  r]°  =  0,  the 


second  term  on  the  left  hand  side  of  (4.4)  becomes 

£  a M"+i(ur, 

n= 0 

=  -k  E  At 

^  n=0 

=  E  jf;  -  £  am-”+1k”, 

+  ^  {  E  A((A-"(«r>, ,”,,”)  -  A-”+‘K”, 

EE  A-'vK'v-‘,,A',,N)-'£1A(f”, 

n=0 

where 

t r  =  («•(•, (",<”-■  )-«•(., (”+i,ur))v,-.v,"rfx,  i< 

r4°  =  o. 

Combining  (4.4)  and  (4.8)  we  see  that 

£  ^  IIMlin  +  ^(uf-1,^,^) 

n=0 

=  EA( (Eh” ■ 

n=0  \A:=1  ' 

We  now  estimate  the  right  hand  side  terms  of  (4.10)  and  for 


(4.8) 

n  <  N  —  1 
(4.9) 


(4.10) 
ease  of 
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notation  adopt  again  the  norm  |||  ■  |||,  where 

|W|||J  =  £  Af  ||<rllln  •  (4.11) 

n=0 

By  the  Schwartz  inequality  we  have 

y>((«\^")<iu«i!i-  nielli  (4.i2) 

n=0 

and 

t(E  '  ,6tij  )  <  At  ||utt||£;2((0i7);£i2(o)  •  (4-13) 

n=0 

Now  using  summation  by  parts  we  have  that 


N- 1 

E  A tT?  = 

n= 0 

=  XA<  J. ’  (o'(x, -o*(x,(“+1,«"))Vfirfl -V(f, 9") 

71  =  0 

=  E{-/  A^([a*(x,r,<"-1)-a*(x,r,u"-1)]V<)-V77"dx} 

n= o  **  -/n  J 

+  l  (a^x,^,^-1)  -a,(x,^,uJV-1))  Vu^-  Vr/^dx 

=  E  -  [jU  (a-(x,iB,ur1)  -  <T(x,<n,un-1))  V«X  ■  V^dx 

+  jf  (a*(x,«n,<n“1)-a*(x,r,un-1))  V^-V^dx  At 
+  J  (a*(x,  tN,  u^-1)  -  a*(x,  tN,  /-*))  Vuf  ■  V^dx  .  (4.14) 


We  note  that 
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=  (^(x,(”,p”-1)K”-'-u”-1)j 

i  C/U  in-fl  -n\c  /„*n— 1  *  .n— 1  \ 

+  )MuA  -  u  ) 

for  any  vn~x  lying  between  u™~1  and  un_1. 

Using  the  Lipschitz  continuity  of  a*  and  the  Schwarz  inequality  we  deduce 
that 

N- 1 

E  AfT? 

n=0 

<  E4‘f  4«lr  HV*fc|U„  (Ilx”l!„,„  +  IICIU) 

n= 1  \  UU  oo,Q 

+  |  llvfi4|L,„(||a„'-'||o!i  +  ||a1r-1||on) 

OO 

+  |  iivw;iL,„(||x"-10in  +  ||r-,||J)iiv,”ii0,1) 

OO  ' 

+  (  ^_1  +  VN_1  )  Vu^ll  Iv^ll 

Vs  o,n  1  o,n /  h  lloo.o  II  '  llo,o 

<  c1(|||,|||  +  |||{|||  +  nielli  +  |||«|||)-HIVxlll 

where  Cx  depends  on  ^  ^  §j|  ,  and  ||V<5tu£||. 
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We  have 


N- 1 

E  A^n  = 

n=0 

N- 1  r 

=  E  A*  /  (<**(*, *n+1X)  -  «(x> *n+\  un+1))  v<+1  •  vdtf 

=  E  At  (jf  (a*(x,tn+\un)  -  a*(x,tn+1un+1))  Vu£+1  •  V^t/ 
+  J  (a*(x,fn+1,un+1)  -  a(x,  <n+1,  un+1))  Vunh+1  ■  V<9,7/n) 

Again  we  apply  summation  by  parts  which  yields 
N- 1 

E  Atf2n  = 

71=0 

(  r  3(1*  \ 

=  £  -(At)2  [JnS>(^-  W-'VuZj  ■  V,Vx 

-  Af  /  ^(-  ,<'v,«'v-1)61u'v-'VS^  •  VijNdx 

Jn  ou 
N- 1 

-  E 

71=1 

+  ^(a*  -  a)(-  ,  t‘ ^/)Vuf  •  Vr/Ndx  . 

Thus,  from  (4.17) 


/  5,  ((a*  -  a)nV<)  •  VVndx 
J  Q 


.(4.16) 


(4.17) 
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N-l 

£  A^2n  < 


n=0 


C  At 


da * 


du 

+At 

II 

+ 

-f-  Ai 

+ 

+  ||o*-fl| 


IVfijtl 


I  du 2 


oc.n  dt2 


d2a' 


dtdu 


du 


dt 


L»((0,T);L»(n)) 

•  l|v^lL,n 


da* 

du 

oo 

da* 

du 

oo 

y- 

-a) 

fi.tt”-1! 


1 


L*((0,T);L*(n)) 

l|V(«O.IL,n |  ll|V»||| 


L2((0,T),L2(fi)) 


dt 


0,0 


||VflN||  ||Vu»||t 


^lia'vsk 

L2((0 ,T);L2(0))  €=0 


IIIVi/lll 


L°°((0,T);L2(O)) 


l|Vi»IL  ■  ||v," 


(4.18) 


We  have 

N-l 

E  A'h”  < 

71  =  0 

<  hi/ -  mi  iimii 

<  c(llM!l  +  M  +  />,+p  1  IIuIIl2((o,t),h?(o)))  ■  • 

(4.19) 
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In  addition  we  see  that 


N-l 

£  A tT?  < 

n=0 

<  "f  Ai(||WIU  +  6-£-  I  A'2||’?“llo.nllv’f"llo,!i  ■  (4.20) 

n=0  \  01  0,0/ 

Applying  results  from  Theorem  3.1  for  q  =  2  and  p  =  1  we  have  that 

Mo,n  <C(h2  +  At). 

Letting  At  =  0(h2)  we  have  that 
N- 1 

£  A tf4n  < 

n=0 

<  (llMII  +  § 

\  0,1  L2((0,T);L2(n)) 

+  A2  ^  )  IIIVolll  .  (4.21) 

Ul  L2(0,T),//2(n))/ 

Combining  the  estimates  (4.15),  (4.18),  (4.19)  and  (4.21)  we  conclude 
that 
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Ill^lll  +  A'N(uhN-\VN,VN) 

—  Ill^^lll  lll^t^ll I  +  At  !|u((||£2((0 ,t);LHQ))  III^ij?III 

+  C,  (III,  III  +  l!lf]||  +  infill  +  I||«|||)I||V,  III 

+  c-(lk'v-,IL+ll"A'-,UlKL 

+  C^At  ||^||//2((0,r);L2(n))  lll^lll 

+  ^2  (lIMII  +  llslll  +  h2  ||wj|L2((0tT);H2(Q))) 

+  IIMII+  § 

Ul  L2((0  ,T);L2(n)) 

+  />!  ^  |||V,||| 

°l  £,2((0,X);//2(Q)) 

where  Cx  and  C2  depend  on  || |  ^  n  and  ||  Vu^l^.  C2  also 

depends  on  UVu/JI^. 

Using  approximation  theory  and  nonlinear  results  from  Theorem  3.1  and 
Lemma  3.1,  we  have 

I*?"!  Tlll^HI  <C(h2  +  At).  (4.22) 

Thus  for  At  —  0(h 2)  we  have 
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Ill^lll  +  A'(uf-\VN,r1N) 

<  e  HI^III2  +  C(e)\h4  +  (At)2  +  6  Vr]N  2|  , 

. 

which  implies  that 

ML  <c(h2+At). 

1,00,  U 

This  provides  the  Too-time  estimate  for  the  backward  difference  procedure. 
As  stated  earlier,  the  generalization  to  the  case  0  <  0  <  |  follows  similar 
arguments.  We  finally  state  these  results  as  a  theorem. 

Theorem  4.1  Let  u(x,t)  be  the  solution  of  problem  (3-4)  and  u£n(x)  be  the 
solution  of  (3.8)  with  0  =  0,  where  *  denotes  numerical  integration  using 
either  interpolation  of  coefficients  or  Gauss  quadrature.  It  is  assumed  that 
q  =  2  and  p  —  1  in  (1.12).  Then 

||Vu(-  ,fn)  —  Vu^R"jon  <  C(h2  +  At0) 

for  any  n  <  N  where  (3  =  1  if  0  <  0  <  \  and  /3  =  2  if  6  =  | . 
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